\(\int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx\) [644]

Optimal result

Integrand size = 23, antiderivative size = 227 \[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=-\frac {2 a b \operatorname {Hypergeometric2F1}\left (2,\frac {m}{2},\frac {2+m}{2},\frac {b^2 \sec ^2(e+f x)}{a^2+b^2}\right ) (d \sec (e+f x))^m}{\left (a^2+b^2\right )^2 f m}+\frac {\operatorname {AppellF1}\left (\frac {1}{2},2,1-\frac {m}{2},\frac {3}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right ) (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2} \tan (e+f x)}{a^2 f}+\frac {b^2 \operatorname {AppellF1}\left (\frac {3}{2},2,1-\frac {m}{2},\frac {5}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right ) (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2} \tan ^3(e+f x)}{3 a^4 f} \]

[Out]

-2*a*b*hypergeom([2, 1/2*m],[1+1/2*m],b^2*sec(f*x+e)^2/(a^2+b^2))*(d*sec(f*x+e))^m/(a^2+b^2)^2/f/m+AppellF1(1/
2,2,1-1/2*m,3/2,b^2*tan(f*x+e)^2/a^2,-tan(f*x+e)^2)*(d*sec(f*x+e))^m*tan(f*x+e)/a^2/f/((sec(f*x+e)^2)^(1/2*m))
+1/3*b^2*AppellF1(3/2,2,1-1/2*m,5/2,b^2*tan(f*x+e)^2/a^2,-tan(f*x+e)^2)*(d*sec(f*x+e))^m*tan(f*x+e)^3/a^4/f/((
sec(f*x+e)^2)^(1/2*m))

Rubi [A] (verified)

Time = 0.31 (sec) , antiderivative size = 227, normalized size of antiderivative = 1.00, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.261, Rules used = {3593, 771, 440, 455, 70, 524} \[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\frac {\tan (e+f x) \sec ^2(e+f x)^{-m/2} (d \sec (e+f x))^m \operatorname {AppellF1}\left (\frac {1}{2},2,1-\frac {m}{2},\frac {3}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right )}{a^2 f}-\frac {2 a b (d \sec (e+f x))^m \operatorname {Hypergeometric2F1}\left (2,\frac {m}{2},\frac {m+2}{2},\frac {b^2 \sec ^2(e+f x)}{a^2+b^2}\right )}{f m \left (a^2+b^2\right )^2}+\frac {b^2 \tan ^3(e+f x) \sec ^2(e+f x)^{-m/2} (d \sec (e+f x))^m \operatorname {AppellF1}\left (\frac {3}{2},2,1-\frac {m}{2},\frac {5}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right )}{3 a^4 f} \]

[In]

Int[(d*Sec[e + f*x])^m/(a + b*Tan[e + f*x])^2,x]

[Out]

(-2*a*b*Hypergeometric2F1[2, m/2, (2 + m)/2, (b^2*Sec[e + f*x]^2)/(a^2 + b^2)]*(d*Sec[e + f*x])^m)/((a^2 + b^2
)^2*f*m) + (AppellF1[1/2, 2, 1 - m/2, 3/2, (b^2*Tan[e + f*x]^2)/a^2, -Tan[e + f*x]^2]*(d*Sec[e + f*x])^m*Tan[e
+ f*x])/(a^2*f*(Sec[e + f*x]^2)^(m/2)) + (b^2*AppellF1[3/2, 2, 1 - m/2, 5/2, (b^2*Tan[e + f*x]^2)/a^2, -Tan[e
+ f*x]^2]*(d*Sec[e + f*x])^m*Tan[e + f*x]^3)/(3*a^4*f*(Sec[e + f*x]^2)^(m/2))

Rule 70

Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(b*c - a*d)^n*((a + b*x)^(m + 1)/(b^(
n + 1)*(m + 1)))*Hypergeometric2F1[-n, m + 1, m + 2, (-d)*((a + b*x)/(b*c - a*d))], x] /; FreeQ[{a, b, c, d, m
}, x] && NeQ[b*c - a*d, 0] &&  !IntegerQ[m] && IntegerQ[n]

Rule 440

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[a^p*c^q*x*AppellF1[1/n, -p,
-q, 1 + 1/n, (-b)*(x^n/a), (-d)*(x^n/c)], x] /; FreeQ[{a, b, c, d, n, p, q}, x] && NeQ[b*c - a*d, 0] && NeQ[n
, -1] && (IntegerQ[p] || GtQ[a, 0]) && (IntegerQ[q] || GtQ[c, 0])

Rule 455

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.), x_Symbol] :> Dist[1/n, Subst[Int
[(a + b*x)^p*(c + d*x)^q, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p, q}, x] && NeQ[b*c - a*d, 0] && EqQ[m
- n + 1, 0]

Rule 524

Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[a^p*c^q*
((e*x)^(m + 1)/(e*(m + 1)))*AppellF1[(m + 1)/n, -p, -q, 1 + (m + 1)/n, (-b)*(x^n/a), (-d)*(x^n/c)], x] /; Free
Q[{a, b, c, d, e, m, n, p, q}, x] && NeQ[b*c - a*d, 0] && NeQ[m, -1] && NeQ[m, n - 1] && (IntegerQ[p] || GtQ[a
, 0]) && (IntegerQ[q] || GtQ[c, 0])

Rule 771

Int[((d_) + (e_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Int[ExpandIntegrand[(a + c*x^2)^p, (d/(d
^2 - e^2*x^2) - e*(x/(d^2 - e^2*x^2)))^(-m), x], x] /; FreeQ[{a, c, d, e, p}, x] && NeQ[c*d^2 + a*e^2, 0] &&
!IntegerQ[p] && ILtQ[m, 0]

Rule 3593

Int[((d_.)*sec[(e_.) + (f_.)*(x_)])^(m_.)*((a_) + (b_.)*tan[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Dist[d^(2*
IntPart[m/2])*((d*Sec[e + f*x])^(2*FracPart[m/2])/(b*f*(Sec[e + f*x]^2)^FracPart[m/2])), Subst[Int[(a + x)^n*(
1 + x^2/b^2)^(m/2 - 1), x], x, b*Tan[e + f*x]], x] /; FreeQ[{a, b, d, e, f, m, n}, x] && NeQ[a^2 + b^2, 0] &&
!IntegerQ[m/2]

Rubi steps \begin{align*} \text {integral}& = \frac {\left ((d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2}\right ) \text {Subst}\left (\int \frac {\left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{(a+x)^2} \, dx,x,b \tan (e+f x)\right )}{b f} \\ & = \frac {\left ((d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2}\right ) \text {Subst}\left (\int \left (\frac {a^2 \left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{\left (a^2-x^2\right )^2}-\frac {2 a x \left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{\left (a^2-x^2\right )^2}+\frac {x^2 \left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{\left (-a^2+x^2\right )^2}\right ) \, dx,x,b \tan (e+f x)\right )}{b f} \\ & = \frac {\left ((d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2}\right ) \text {Subst}\left (\int \frac {x^2 \left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{\left (-a^2+x^2\right )^2} \, dx,x,b \tan (e+f x)\right )}{b f}-\frac {\left (2 a (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2}\right ) \text {Subst}\left (\int \frac {x \left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{\left (a^2-x^2\right )^2} \, dx,x,b \tan (e+f x)\right )}{b f}+\frac {\left (a^2 (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2}\right ) \text {Subst}\left (\int \frac {\left (1+\frac {x^2}{b^2}\right )^{-1+\frac {m}{2}}}{\left (a^2-x^2\right )^2} \, dx,x,b \tan (e+f x)\right )}{b f} \\ & = \frac {\operatorname {AppellF1}\left (\frac {1}{2},2,1-\frac {m}{2},\frac {3}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right ) (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2} \tan (e+f x)}{a^2 f}+\frac {b^2 \operatorname {AppellF1}\left (\frac {3}{2},2,1-\frac {m}{2},\frac {5}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right ) (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2} \tan ^3(e+f x)}{3 a^4 f}-\frac {\left (a (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2}\right ) \text {Subst}\left (\int \frac {\left (1+\frac {x}{b^2}\right )^{-1+\frac {m}{2}}}{\left (a^2-x\right )^2} \, dx,x,b^2 \tan ^2(e+f x)\right )}{b f} \\ & = -\frac {2 a b \operatorname {Hypergeometric2F1}\left (2,\frac {m}{2},\frac {2+m}{2},\frac {b^2 \sec ^2(e+f x)}{a^2+b^2}\right ) (d \sec (e+f x))^m}{\left (a^2+b^2\right )^2 f m}+\frac {\operatorname {AppellF1}\left (\frac {1}{2},2,1-\frac {m}{2},\frac {3}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right ) (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2} \tan (e+f x)}{a^2 f}+\frac {b^2 \operatorname {AppellF1}\left (\frac {3}{2},2,1-\frac {m}{2},\frac {5}{2},\frac {b^2 \tan ^2(e+f x)}{a^2},-\tan ^2(e+f x)\right ) (d \sec (e+f x))^m \sec ^2(e+f x)^{-m/2} \tan ^3(e+f x)}{3 a^4 f} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains complex when optimal does not.

Time = 18.79 (sec) , antiderivative size = 2453, normalized size of antiderivative = 10.81 \[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\text {Result too large to show} \]

[In]

Integrate[(d*Sec[e + f*x])^m/(a + b*Tan[e + f*x])^2,x]

[Out]

((d*Sec[e + f*x])^m*((-2*a*b*(-1 + (Sec[e + f*x]^2)^(m/2)))/m + (a^2 - b^2)*Hypergeometric2F1[1/2, 1 - m/2, 3/
2, -Tan[e + f*x]^2]*Tan[e + f*x] + (2*a*b*AppellF1[-m, -1/2*m, -1/2*m, 1 - m, (a - I*b)/(a + b*Tan[e + f*x]),
(a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2))/(m*((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2
)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)) + (b*(a^2 + b^2)*AppellF1[1 - m, -1/2*m, -1/2*m, 2 - m,
(a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2))/((-1 + m)*((b*(-I + T
an[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*(a + b*Tan[e + f
*x]))))/(f*(a + b*Tan[e + f*x])^2*((a^2 - b^2)*Hypergeometric2F1[1/2, 1 - m/2, 3/2, -Tan[e + f*x]^2]*Sec[e + f
*x]^2 - 2*a*b*(Sec[e + f*x]^2)^(m/2)*Tan[e + f*x] + (2*a*b*AppellF1[-m, -1/2*m, -1/2*m, 1 - m, (a - I*b)/(a +
b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2)*Tan[e + f*x])/(((b*(-I + Tan[e + f*x])
)/(a + b*Tan[e + f*x]))^(m/2)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)) - (b^2*(a^2 + b^2)*AppellF1
[1 - m, -1/2*m, -1/2*m, 2 - m, (a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2
)^(1 + m/2))/((-1 + m)*((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*((b*(I + Tan[e + f*x]))/(a + b*Tan
[e + f*x]))^(m/2)*(a + b*Tan[e + f*x])^2) + (b*(a^2 + b^2)*m*AppellF1[1 - m, -1/2*m, -1/2*m, 2 - m, (a - I*b)/
(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2)*Tan[e + f*x])/((-1 + m)*((b*(-I +
Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*(a + b*Tan[e +
f*x])) + (2*a*b*(Sec[e + f*x]^2)^(m/2)*(-1/2*((a - I*b)*b*m^2*AppellF1[1 - m, 1 - m/2, -1/2*m, 2 - m, (a - I*
b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*Sec[e + f*x]^2)/((1 - m)*(a + b*Tan[e + f*x])^2) - ((
a + I*b)*b*m^2*AppellF1[1 - m, -1/2*m, 1 - m/2, 2 - m, (a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e
+ f*x])]*Sec[e + f*x]^2)/(2*(1 - m)*(a + b*Tan[e + f*x])^2)))/(m*((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x])
)^(m/2)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)) + (b*(a^2 + b^2)*(Sec[e + f*x]^2)^(m/2)*(((a - I*
b)*b*(1 - m)*m*AppellF1[2 - m, 1 - m/2, -1/2*m, 3 - m, (a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e
+ f*x])]*Sec[e + f*x]^2)/(2*(2 - m)*(a + b*Tan[e + f*x])^2) + ((a + I*b)*b*(1 - m)*m*AppellF1[2 - m, -1/2*m, 1
- m/2, 3 - m, (a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*Sec[e + f*x]^2)/(2*(2 - m)*(a +
b*Tan[e + f*x])^2)))/((-1 + m)*((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*((b*(I + Tan[e + f*x]))/(
a + b*Tan[e + f*x]))^(m/2)*(a + b*Tan[e + f*x])) - (a*b*AppellF1[-m, -1/2*m, -1/2*m, 1 - m, (a - I*b)/(a + b*T
an[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2)*((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f
*x]))^(-1 - m/2)*(-((b^2*Sec[e + f*x]^2*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x])^2) + (b*Sec[e + f*x]^2)/(a +
b*Tan[e + f*x])))/((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2) - (b*(a^2 + b^2)*m*AppellF1[1 - m, -1/2
*m, -1/2*m, 2 - m, (a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2)*((b*
(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(-1 - m/2)*(-((b^2*Sec[e + f*x]^2*(-I + Tan[e + f*x]))/(a + b*Tan[e
+ f*x])^2) + (b*Sec[e + f*x]^2)/(a + b*Tan[e + f*x])))/(2*(-1 + m)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x
]))^(m/2)*(a + b*Tan[e + f*x])) - (a*b*AppellF1[-m, -1/2*m, -1/2*m, 1 - m, (a - I*b)/(a + b*Tan[e + f*x]), (a
+ I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2)*((b*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(-1 - m/2)*(
-((b^2*Sec[e + f*x]^2*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x])^2) + (b*Sec[e + f*x]^2)/(a + b*Tan[e + f*x])))/
((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2) - (b*(a^2 + b^2)*m*AppellF1[1 - m, -1/2*m, -1/2*m, 2 - m,
(a - I*b)/(a + b*Tan[e + f*x]), (a + I*b)/(a + b*Tan[e + f*x])]*(Sec[e + f*x]^2)^(m/2)*((b*(I + Tan[e + f*x])
)/(a + b*Tan[e + f*x]))^(-1 - m/2)*(-((b^2*Sec[e + f*x]^2*(I + Tan[e + f*x]))/(a + b*Tan[e + f*x])^2) + (b*Sec
[e + f*x]^2)/(a + b*Tan[e + f*x])))/(2*(-1 + m)*((b*(-I + Tan[e + f*x]))/(a + b*Tan[e + f*x]))^(m/2)*(a + b*Ta
n[e + f*x])) + (a^2 - b^2)*Sec[e + f*x]^2*(-Hypergeometric2F1[1/2, 1 - m/2, 3/2, -Tan[e + f*x]^2] + (1 + Tan[e
+ f*x]^2)^(-1 + m/2))))

Maple [F]

\[\int \frac {\left (d \sec \left (f x +e \right )\right )^{m}}{\left (a +b \tan \left (f x +e \right )\right )^{2}}d x\]

[In]

int((d*sec(f*x+e))^m/(a+b*tan(f*x+e))^2,x)

[Out]

int((d*sec(f*x+e))^m/(a+b*tan(f*x+e))^2,x)

Fricas [F]

\[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\int { \frac {\left (d \sec \left (f x + e\right )\right )^{m}}{{\left (b \tan \left (f x + e\right ) + a\right )}^{2}} \,d x } \]

[In]

integrate((d*sec(f*x+e))^m/(a+b*tan(f*x+e))^2,x, algorithm="fricas")

[Out]

integral((d*sec(f*x + e))^m/(b^2*tan(f*x + e)^2 + 2*a*b*tan(f*x + e) + a^2), x)

Sympy [F]

\[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\int \frac {\left (d \sec {\left (e + f x \right )}\right )^{m}}{\left (a + b \tan {\left (e + f x \right )}\right )^{2}}\, dx \]

[In]

integrate((d*sec(f*x+e))**m/(a+b*tan(f*x+e))**2,x)

[Out]

Integral((d*sec(e + f*x))**m/(a + b*tan(e + f*x))**2, x)

Maxima [F]

\[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\int { \frac {\left (d \sec \left (f x + e\right )\right )^{m}}{{\left (b \tan \left (f x + e\right ) + a\right )}^{2}} \,d x } \]

[In]

integrate((d*sec(f*x+e))^m/(a+b*tan(f*x+e))^2,x, algorithm="maxima")

[Out]

integrate((d*sec(f*x + e))^m/(b*tan(f*x + e) + a)^2, x)

Giac [F]

\[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\int { \frac {\left (d \sec \left (f x + e\right )\right )^{m}}{{\left (b \tan \left (f x + e\right ) + a\right )}^{2}} \,d x } \]

[In]

integrate((d*sec(f*x+e))^m/(a+b*tan(f*x+e))^2,x, algorithm="giac")

[Out]

integrate((d*sec(f*x + e))^m/(b*tan(f*x + e) + a)^2, x)

Mupad [F(-1)]

Timed out. \[ \int \frac {(d \sec (e+f x))^m}{(a+b \tan (e+f x))^2} \, dx=\int \frac {{\left (\frac {d}{\cos \left (e+f\,x\right )}\right )}^m}{{\left (a+b\,\mathrm {tan}\left (e+f\,x\right )\right )}^2} \,d x \]

[In]

int((d/cos(e + f*x))^m/(a + b*tan(e + f*x))^2,x)

[Out]

int((d/cos(e + f*x))^m/(a + b*tan(e + f*x))^2, x)